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We report the first observation of diffractively produced W bosons. In a sample 
of W — * ev events produced in pp collisions at JH=1.8 TeV, we find an excess 
of events with a forward rapidity gap, which is attributed to diffraction. The 
probability that this excess is consistent with non-diffractive production is 1.1 x 
10~ 4 (3.8a). The relatively low fraction of W + Jet events observed within this 
excess implies that mainly quarks from the pomeron, which mediates diffraction, 
participate in W production. The diffractive to non-diffractive W production ratio 
is found to be R w = (1.15 ± 0.55)%. 

O 

Approximately 15% of high energy pp inelastic collisions are due to single diffraction 

O 



dissociation, a process in which the incident p or p escapes intact losing a fraction 
■ £ < 0.1 of its initial forward momentum. Experiments have shown that the leading 

role in diffraction is played by the pomeron 0, which carries the quantum numbers of 
the vacuum. In QCD the pomeron is a colorless entity, whose exchange in an event is 
marked by a "rapidity gap", i.e. a large region of pseudorapidity || devoid of particles. 

The partonic structure of the pomeron was first investigated by the UA8 experiment 
[3], which studied diffractive dijet production at the CERN SppS collider at yfs = 630 
GeV, and more recently by the HI |], [F| and ZEUS || |9| experiments in diffractive 
deep inelastic scattering (DDIS) || [?| |8[] and dijet photoproduction |J in ep collisions at 
y/s ~ 300 GeV at HERA. All experiments find that a substantial fraction of the pomeron 
structure is "hard", i.e. consists of partons carrying a large fraction of the pomeron 
momentum. From the DDIS experiments, which probe directly the quark component of 
the pomeron, the hard-quark component is estimated to account for about one third of 
the pomeron momentum. At the Tevatron pp collider, a hard-quark pomeron structure 
would lead to detectable diffractive W production [|10j, which to leading order occurs 
through q'q — > W. For a hard-gluon dominated pomeron, W production can occur 
through qg — > Wq', but at a rate lower by order a s and always in association with a jet. 

In this paper, we present the results of a measurement of diffractive W production 
in pp collisions at \/s = 1.8 TeV using the CDF detector at the Tevatron. Diffraction is 
tagged by the presence of a rapidity gap in an event in association with the following ex- 
pected characteristic features. In a diffractive W — > e ± i > event produced in a p collision 
with a pomeron (V) emitted by the proton, the rapidity gap is expected to be at posi- 
tive r] (p-direction) and the lepton boosted towards negative r] (angle-gap correlation). 
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Also, since the pomeron is quark-flavor symmetric, and since from energy considerations 
mainly valence quarks from the p participate in producing the W, approximately twice 
as many electrons as positrons are expected (charge-gap correlation). These correlations 
can be seen in the Monte Carlo (MC) generated distributions of Fig. 1. The opposite 
correlations are, of course, expected for p — V collisions with the pomeron emitted by 
the p. In non-diffractive events, where rapidity gaps may arise from fluctuations in the 
event particle multiplicity, MC simulations using the PYTHIA |TT| program show that 
there are no significant angle-gap or charge-gap correlations. 

We simulate diffractive events using the POMPYT |Oj MC program, which is based 



on the Ingelman-Schlein model for hard diffraction [13|]. The cross section for pp — > pX 
may be written as 

where K is a constant, £ is the fraction of the momentum of the proton carried by the 
pomeron, t is the square of the four- momentum transfer, a(t) = l+e+a't is the pomeron 
trajectory, F(t) the nucleon form factor, s = £s the center of mass energy squared in 
the V — p reference frame, and 0j^(s) the V — p total cross section. This equation 
suggests the interpretation of single diffraction dissociation as a process in which a flux of 
pomerons, /p/ p (£,i), emitted by the proton interacts with the antiproton. This concept 
of factorization was extended Jl3| to hard processes by treating the "pomeron flux factor" 
as a flux of particle-like pomerons with a unique partonic structure. In POMPYT, the 
collision of this flux of pomerons with the nucleon is handled by PYTHIA. All our MC 
simulations include a simulation of the CDF detector. 



The CDF detector is described in detail elsewhere [|14], [T5|]. In the rapidity gap 
analysis we use the "beam-beam counters" (BBC) and the forward electromagnetic 
(EM) and hadronic (HA) calorimeters. The BBC [14| consist of a square array of 16 
scintillation counters on each side of the interaction point covering approximately the 
region 3.2 < \rj\ < 5.9. The forward calorimeters cover the region 2.4 < \t]\ < 4.2 and 
have projective tower geometry with tower size At] ■ A0 = 0.1 x 5°, where is the 
azimuthal angle. An energy threshold of 1.5 GeV (sum of EM plus HA energies) is used 
for each tower to exclude calorimeter noise. 

The data sample was obtained during collider runs 1A (1992-1993) and IB (1994- 
1995) by triggering on an electron of high transverse momentum, Pt = P sin# e , and on 
missing transverse energy, Et [0. We used events with Et > 20 GeV and an isolated 
ji"7| electron of Et > 20 GeV in the central region, \rj\ < 1.1, where the tracks of charged 



particles can be completely reconstructed. After implementing a cut retaining events 
with one primary vertex only, 8246 events remained. The one-vertex cut was imposed to 
exclude events with two interactions in the same beam-beam crossing, since the overlay 
of a "minimum bias" on a diffractive W event could eliminate the rapidity gap. 

We search for a diffractive W signal by analyzing the correlations between the r\ 
of the electron, rj e , or the sign of its charge, C e , and the multiplicity of one or the 
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other of the BBCs. Each event enters into two distributions, one with r/ e ■ t/bbc < 
(angle-correlated) or C e ■ t/bbc < (charge-correlated), and the other with r\ e ■ t/bbc > 
(angle-anticorrelated) or C e ■ t/bbc > (charge-anticorrelated). A doubly-correlated 
(anticorrelated) distribution is the BBC multiplicity distribution for events with T/ e -C e > 
and r/ e ■ t/bbc < (r) e ■ t/bbc > 0). Fig. 2 shows the observed correlations as a function 
of BBC multiplicity, Nbbc, for events with tower multiplicity, Nt, less than 8 in the 
forward calorimeter adjacent to a given BBC. The cut on Nt is imposed to reduce the 
non-diffractive contribution to the signal, since the signal is concentrated at low N B bc 
and is expected to have low N T as well. Fig. 2a shows the angle and charge doubly- 
correlated (solid) and doubly-anticorrelated (dashed) BBC multiplicities. The peaking 
at high multiplicities is caused by saturation due to the finite BBC segmentation. The 
two distributions agree well above the first three bins, but the correlated distribution 
has an excess in the first bin, consistent with the signature expected from diffractive 
events with a rapidity gap. This excess can be seen more clearly in Fig. 2b, which 
shows the bin-by-bin asymmetry (difference divided by sum) of the two distributions of 
Fig. 2a. An excess is also seen in the individual angle (Fig. 2c) and charge (Fig. 2d) 
correlated asymmetries, as expected for diffractive production. From MC simulations 
of non-diffractive W production and using Poisson statistics, the probability that the 
observed excess in the first bin of both the angle and charge correlated distributions is 
due to simultaneous fluctuations in the non-diffractive background was estimated to be 
1.1 x 10" 4 . 

The quark to gluon fraction of the partons of the pomeron participating in W pro- 
duction may be evaluated from the fraction of diffractive W + Jet events observed. 
Simulations performed with a hard-gluon (quark) pomeron structure predict the frac- 
tion of diffractive W events containing at least one jet with Et > 6 GeV (within an rj — <fi 
cone radius of 0.7) to be 0.66 (0.20). For non-diffractive W events with similar kinemat- 
ics the predicted "jet fraction" is 0.34, consistent with measurements in a non-diffractive 
data sample. In the first bin of Fig. 2a (solid histogram) there are 34 events, among 
which we estimate 21 to be diffractive and 13 non-diffractive. Multiplying these numbers 
by the corresponding predicted jet fractions yields an expectation of 18.4 ±2.8 (8.8 ±2.5) 
events with a jet for a hard-gluon (quark) pomeron structure. The data contain 8 events 
with a jet, which implies that predominantly quarks from the pomeron participate in W 
production. 

We use the doubly-correlated distributions of Fig. 2a to evaluate the ratio, R, of 
diffractive to non-diffractive W production rates. As a ratio, R is insensitive to lepton 
selection cuts or inefficiencies and to the uncertainty in the luminosity. The acceptance 
for diffractive events is obtained from POMPYT using a hard-quark pomeron structure 
of the form (3G(f3) = 6(3(1 — (3), where /3 is the fraction of the momentum of the pomeron 
carried by the quark. In order to check for possible systematic effects due to BBC noise or 
inefficiencies that could distort the low multiplicity binning and thereby give an incorrect 
R ratio, we evaluate R using events with a BBC multiplicity upper bound, Nb, and we 
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vary Nb from zero to seven. Fig. 3a shows the resulting R values, and Fig. 3b the MC 
"gap-acceptance", as a function of Nb- The gap-acceptance for bin Nb is defined as the 
fraction of events with N B bc < N B (the lepton acceptance is not included here). The 
errors in the points of Fig. 3a, which are statistical, increase with increasing Nb as more 
background is being subtracted. To reduce the sensitivity of the result to the acceptance 
calculation, we retain as our signal the value R = (1.03 ± 0.46)% of the Nb = 2 bin, 
where the acceptance is 81% and varies relatively slowly with Nb- 

As a systematic uncertainty in the gap-acceptance calculation we assign ±13%, which 
is one half of the difference between the acceptances of Nb = 1 and Nb = 3 divided 
by the acceptance of Nb = 2. In deriving the ratio R we assumed that the non- 
diffractive contributions to the correlated and anticorrelated distributions in Fig. 2a are 
identical. This assumption is justified by the excellent matching of the two distributions 
for Nb > 3. A possible mismatch of the distributions within the available statistics 
introduces a systematic uncertainty, which was evaluated as follows. We made a straight 
line fit to the asymmetry of bins 4-10 of Fig. 2b, and extrapolated the fit into bins 1-3. 
For each of the bins 1-3, we multiplied the extrapolated asymmetry and its error by 
twice the number of anticorrelated events, since the average number of non-diffractive 
correlated and anticorrelated events is expected to be the same, and added up the results 
for the three bins. Treating the sum as a signal yields a diffractive to non-diffractive 
ratio of (0.01 ± 0.11)%, which is consistent with zero. We treat the error of ±0.11% as 
a systematic uncertainty in our measured value of R and add it in quadrature to the 
gap-acceptance uncertainty to obtain a combined systematic uncertainty of ±0.18%. 

From a study of the rate of W events versus instantaneous luminosity we estimate 
that a correction of 0.95 ± 0.05(syst) must be applied to R due to the different effi- 
ciency of the single vertex cut for diffractive and non-diffractive events. In addition, 
we apply a correction for the BBC occupancy by particles from a second interaction 
that does not have a reconstructed vertex. From a study of a sample of 98000 events 
recorded by triggering the detector on beam-beam crossings only, we determined that 
the probability of finding more than two hits in a BBC is 15%, corresponding to a 
BBC livetime acceptance of 0.85 by which we divide R. The corrected value for R is 
Rw = [1.15 ± 0.51(stat) ± 0.20(syst)]%. From MC simulations we estimate that the 
diffractive events are concentrated at ^-values in the range 0.01-0.05. 

Below we compare our results with POMPYT predictions and with results from other 
experiments. The predictions depend on the assumed pomeron structure function and 
on the form and normalization of the pomeron flux factor, f-p/ p (l;,t). We first use the 
"standard" flux factor |18[ with parameters a(t) = 1.115 + 0.26 1 and K = 0.73 GeV~ 2 ; 



for the nucleon form factor we use @ F(t) = (Am 2 p - 2M)(Am 2 p - fp 1 [1 - t/0.7] . For 
a two (three) flavor hard-quark pomeron structure of the form (3G(f3) = 6/3(1 — (3) we 
obtain i?(y=24% (16%), while for a hard-gluon structure of the same form, R^ = 1.1%. 
Our measured ratio, Rw = (1.15 ± 0.55)%, favors a purely gluonic pomeron, which 
however is incompatible with the low fraction of diffractive W + Jet events we observe. 
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The HERA experiments on DDIS |6|, § at 8.5 < Q 2 < 65 GeV 2 report a quark component 
in the pomeron structure which is flat in /3, rises slowly with Q 2 at any given fixed (3, and 
accounts for a fraction of about one third of the momentum of the pomeron, assuming the 
standard pomeron flux. Independent of the pomeron flux normalization, by combining 
diffractive dijet photoproduction and DDIS results, the ZEUS collaboration reports an 
integrated hard-quark momentum fraction of 0.2 < f q < 0.7, while the HI collaboration 
0, from a QCD analysis of DDIS, obtains /, w 0.2 at Q 2 ~ 60 GeV 2 . The Q 2 evolution 
from Q 2 = 60 GeV 2 to Q 2 = M w of the pomeron structure function proposed by 
HI does not change significantly the quark component participating in W production. 
Using a pomeron with a hard-quark fraction of 0.2 and a gluon fraction of 0.8, POMPYT 
predicts ratios Rw of 5.7% (4.1%) for two (three) quark flavors, which are larger than 
our measured value of (1.15 ± 0.55)% by more than eight (five) standard deviations. 

We now compare our results with POMPYT predictions using the "renormalized" 
pomeron flux [18], defined as the standard flux normalized, if its integral exceeds unity, 
to one pomeron per nucleon. The normalization factor is ~ 9 at y/s = 1.8 TeV (CDF) 
and « 1 at HERA (see [0]). The predictions for R w become 2.7% (1.8%) for a two 
(three) flavor pure hard-quark and 0.12% for a pure hard-gluon pomeron structure. 
Based on these predictions, our R w value of (1.15 ±0.55)% implies hard-quark fractions 
of f q = 0.4 ± 0.2 (0.6 ± 0.3) for two (three) quark flavors. These fractions are consistent 
with the ZEUS and HI results of 0.2 < f q < 0.7 and f q ~ 0.2, respectively. Assuming 
the momentum sum rule, f q + f g = 1, the predicted fractional gluon contribution to 
R w is (0.12%)(1 - / 9 )/[(0.12%)(l - /,) + Af q ], where ,4=2.7% (1.8%) for two (three) 
quark flavors. From our values of f q , the gluon contribution to Rw is predicted to be 
6.6% (4.2%) for two (three) quark flavors, which can explain the low fraction of W + Jet 
events we observe. 

In conclusion, we have observed diffractive W production in pp collisions at ^/s = 1.8 
TeV and measured the ratio of diffractive to non-diffractive production rates to be 
Rw = (1.15 ±0.55)%. The relatively small fraction of diffractive W + Jet events 
we observe implies that mainly quarks from the pomeron participate in diffractive W 
production. 
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Figure 1: Monte Carlo generated 77-distributions: (a) Particle densities for non- 
diffractive (solid) and for diffractive (dashed) W events for pomerons of beam momentum 
fraction £ = 0.03 emitted by protons (at positive 77); the small bump at 77 pa 8.5 is caused 
by the leading protons, (b) Electrons and positrons from diffractive W ± (— > e ± u) events 
for all pomerons of £ < 0.1 emitted by protons (the vertical dashed lines define the 
boundaries of the region of this measurement). 
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Figure 2: (a) Electron angle and charge doubly-correlated (solid) and anticorrelated 
(dashed) distributions (see text) versus BBC multiplicity, and (b) the corresponding 
asymmetry, defined as the bin- by-bin difference over sum of the two distributions in (a). 
The diffractive signal is seen in the first bin as an excess of events in the correlated 
distribution in (a), and as a positive asymmetry in (ft). An asymmetry is also seen in 
the first bin of the individual angle (c) and charge (d) distributions. 
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Figure 3: (a) Diffractive to non-diffractive W production ratio (not corrected for BBC 
occupancy or one-vertex cut efficiency) as a function of upper bound BBC multiplicity, 
N B - The solid line passes through the N B = 2 point, which we use as our result; (b) 
gap-acceptance for angle-gap and charge-gap doubly-correlated (solid) and anticorrelated 
(dashed) diffractive events with an electron within \rj\ < 1.1. 
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